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Abstract The eigenvalues of a matrix polynomial can be determined classically by 
solving a generalized eigenproblem for a linearized matrix pencil, for instance by 
writing the matrix polynomial in companion form. We introduce a general scaling 
technique, based on tropical algebra, which applies in particular to this compan- 
ion form. This scaling, which is inspired by an earlier work of Akian, Bapat, and 
Gaubert, relies on the computation of "tropical roots". We give explicit bounds, in 
a typical case, indicating that these roots provide accurate estimates of the order of 
magnitude of the different eigenvalues, and we show by experiments that this scal- 
ing improves the accuracy (measured by normwise backward error) of the computa- 
tions, particularly in situations in which the data have various orders of magnitude. 
In the case of quadratic polynomial matrices, we recover in this way a scaling due to 
Fan, Lin, and Van Dooren, which coincides with the tropical scaling when the two 
tropical roots are equal. If not, the eigenvalues generally split in two groups, and the 
tropical method leads to making one specific scaling for each of the groups. 



1 Introduction 

A classical problem is to compute the eigenvalues of a matrix polynomial 
P{X)=Ao+Ail + ---+Adl'' 
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where A/ € C"^",l = 0. . .t/ are given. The eigenvalues are defined as the solutions 
of det(f (A)) = 0. If A is an eigenvalue, the associated right and left eigenvectors x 
and y £ C" are the non-zero solutions of the systems P{X)x = and y*P{X) = 0, 
respectively. A common way to solve this problem, is to convert f into a "linearized" 
matrix pencil 

L(A)=Ax + y, xjeC"'"'"'' 

with the same spectrum as P and solve the eigenproblem for L, by standard nu- 
merical algorithms like the QZ method [16]. If D and D' are invertible diagonal 
matrices, and if a is a non-zero scalar, we may consider equivalently the scaled 
pencilDL(aA)D'. 

The problem of finding the good linearizations and the good scaUngs has received 
a considerable attention. The backward error and conditioning of the matrix pencil 
problem and of its linearizations have been investigated in particular in works of 
Tisseur, Li, Higham, and Mackey, see [17, 11, 12]. 

A scaling on the eigenvalue parameter to improve the normwise backward error 
of a quadratic polynomial matrix was proposed by Fan, Lin, and Van Dooren [8]. 
This scaling only relies on the norms yi I ~ 0, 1,2. In this paper, we intro- 

duce a new family of scalings which also rely on these norms. The degree d is now 
arbitrary. 

These scalings originate from the work of Akian, Bapat, and Gaubert [2, 1], 
in which the entries of the matrices A/ are functions, for instance Puiseux series, 
of a (perturbation) parameter f . The valuations (leading exponents) of the Puiseux 
series representing the different eigenvalues were shown to coincide, under some 
genericity conditions, with the points of non-differentiability of the value function of 
a parametric optimal assignment problem (the tropical eigenvalues), a result which 
can be interpreted in terms of amoebas [13]. Indeed, the definition of the tropical 
eigenvalues in [2, 1] makes sense in any field with valuation. In particular, when the 
coefficients belong to C, we can take the map z ^ log |z| from C to IR U {— 0°} as the 
valuation. Then, the tropical eigenvalues are expected to give, again under some non 
degeneracy conditions, the correct order of magnitude of the different eigenvalues. 

The tropical roots used in the present paper are an approximation of the tropical 
eigenvalues, relying only on the norms 7/ = ||A/||. A better scaling may be achieved 
by considering the tropical eigenvalues, but computing these eigenvalues requires 
0{nd) calls to an optimal assignment algorithm, whereas the tropical roots consid- 
ered here can be computed in 0{d) time, see Remark 3 below for more information. 
We examine such extensions in a further work. 

As an illustration, consider the following quadratic polynomial matrix 

By applying the QZ algorithm on the first companion form of P(A) we get the 
eigenvalues -Inf,- 7.731e-19 , Inf, 3.588e-19, by using the scaling proposed in [8] 
we get -Inf, -3.250e-19, Inf, 3.588e-19. However by using the tropical scaling we 
can find the four eigenvalues properly: - 7.250e-18 ± 9.744e-18i, - 2.102eH-17 ± 
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7.387e+17i. The result was shown to be correct (actually, up to a 14 digits pre- 
cision) with PARI, in which an arbitrarily large precision can be set. The above 
computations were performed in Matlab (version 7.3.0). 

The paper is organized as follows. In Section 2, we recall some classical facts of 
max-plus or tropical algebra, and show that the tropical roots of a tropical polyno- 
mial can be computed in linear time, using a convex hull algorithm. Section 3 states 
preliminary results concerning matrix pencils, linearization and normwise backward 
error 

In Section 4, we describe our scaling method. In Section 5, we give a theorem 
locating the eigenvalues of a quadratic polynomial matrix, which provides some the- 
oretical justification of the method. Finally in Section 6, we present the experimental 
results showing that the tropical scaling can highly reduce the normwise backward 
error of an eigenpair We consider the quadratic case in Section 6.1 and the general 
case in Section 6.2. For the quadratic case, we compare our results with the scaling 
proposed in [8]. 



2 Tropical polynomials 

The max-plus semiring Rmax, is the set IR U {— °°}, equipped with max as addition, 
and the usual addition as multiplication. It is traditional to use the notation © for 
max (so 2 3 = 3), and (K) for + (so 1 (g) 1 = 2). We denote by (D the zero element of 
the semiring, which is such that ©©« = «, here (D = and by 1 the unit element 
of the semiring, which is such that l(E)a = a(g)l = a, here 1=0. We refer the reader 
to [4, 14, 3] for more background. 

A variant of this semiring is the max-times semiring Rmax.x, which is the set of 
nonnegative real numbers IR+, equipped with max as addition, and x as multipli- 
cation. This semiring is isomorphic to Rmax by the map x t-^ logx. So, every notion 
defined over Rmax has an Rmax.x analogue that we shall not redefine explicitly. In the 
sequel, the word "tropical" will refer indifferently to any of these algebraic struc- 
tures. 

Consider a max-plus (formal) polynomial of degree n in one variable, i.e., a 
formal expression P = 0o<A:<„^'<.-^*^ in which the coefficients P^. belong to Rmax, 
and the associated numerical polynomial, which, with the notation of the classi- 
cal algebra, can be written as p{x) = maxo<i:<„/'<: + ^x. Cuninghame-Green and 
Meijer showed [7] that the analogue of the fundamental theorem of algebra holds 
in the max-plus setting, i.e., that p{x) can be written uniquely as p{x) = P„ + 
£i<jt<„max(x,cj.), where ci,...,c„ £ Rmax are the roots, i.e., the points at which 
the maximum attained at least twice. This is a special case of more general notions 
which have arisen recently in tropical geometry [13]. The multiplicity of the root 
c is the cardinality of the set {fc e {1, . . . ,n} | q = c}. Define the Newton polygon 
A{P) of P to be the upper boundary of the convex hull of the set of points [k^Pi^), 
k = 0, . . . ,n. This boundary consists of a number of linear segments. An application 
of Legendre-Fenchel duality (see [2, Proposition 2.10]) shows that the opposite of 



4 



Stephane Gaubert and Meisam Sharify 



the slopes of these segments are precisely the tropical roots, and that the multiplicity 
of a root coincides with the horizontal width of the corresponding segment. (Actu- 
ally, min-plus polynomials are considered in [2], but the max-plus case reduces to 
the min-plus case by an obvious change of variable). Since the Graham scan algo- 
rithm [10] allows us to compute the convex hull of a finite set of points by making 
0{n) arithmetical operations and comparisons, provided that the given set of points 
is already sorted by abscissa, we get the following result. 

Proposition 1. The roots of a max-plus polynomial in one variable can be computed 
in linear time. □ 

The case of a max-times polynomial reduces to the max-plus case by replacing 
every coefficient by its logarithm. The exponentials of the roots of the transformed 
polynomial are the roots of the original polynomial. 



3 Matrix pencil and normwise backward error 



Let us come back to the eigenvalue problem for the matrix pencil P{X) = Aq + 

A\X A hA^/A''. There are many ways to construct a "linearized" matrix pencil 

L(A) =XX + Y, X,Y ^ £ndx„d ^jjjj jjjg g^j^g spectrum as P(A), see [15] for a 
general discussion. In particular, the first companion form XXi + Yi is defined by 



Xi =diag(A^, 



Yi 



(Ak-i Ak-2 ■■■ Aq\ 
-I„ ... 



V 



-/„ ; 



In the experimental part of this work, we are using this linearization. 

To estimate the accuracy of a numerical algorithm computing an eigenpair, we 
shall consider, as in [17], the normwise backward error. The latter arises when con- 
sidering a perturbation 

AP^AAo + AAiX + --- + AAdX'' . 

The backward error of an approximate eigenpair (i, A ) of P is defined by 

Tj(i,A) =min{e: {P{l)+AP{l))x = Q,\\AAi\\2 < e||£/||2,/ = 0, . . .m} . 

The matrices £/ representing tolerances. The following computable expression for 
J] (i, A ) is given in the same reference. 



Tj(i,A) 



a||i||2 



where r = P{X)x and a = ^ |A |'|[£'/||2. In the sequel, we shall take £/ = A/ 
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Our aim is to reduce the normwise backward error, by a scaling of the eigenvalue 
A = ajJ., where a is the scaling parameter This kind of scaling for quadratic poly- 
nomial matrix was proposed by Fan, Lin and Van Dooren [8]. We next introduce a 
new scaling, based on the tropical roots. 



4 Construction of the tropical scaling 

Consider the matrix pencil modified by the substitution A = ajJ. 

P(Ai)=Ao+AiAi + ---+Arf/ 

where A, ~ jSa'A,. 

The tropical scaling which we next introduce is characterized by the property 
that a and j3 are such that P{lJ.) has at least two matrices A, with an (induced) 
Euclidean norm equal to one, whereas the Euclidean norm of the other matrices are 
all bounded by one. This scaling is inspired by the work of M. Akian and R. Bapat 
and S. Gaubert [1], which concerns the perturbation of the eigenvalues of a matrix 
pencil. The theorem on the location of the eigenvalues which is stated in the next 
section provides some justification for the present scaling. 

We associate to the original pencU the max-times polynomial 

tp(jc) =max(7o,7iA,--- ,7dA'') , 

where 

Yi--=\\Ai\\ 

(the symbol t stands for "tropical"). Let ai < < ■ ■ ■ < OCii he the tropical roots of 
tp{x) counted with multiplicities. For each a,-, the maximum is attained by at least 
two mononomials. Subsequently, the transformed polynomial q{x) :— I5jtp{aix), 
with j3, :— (tp(a,))^' has two coefficients of modulus one, and all the other co- 
efficients have modulus less than or equal to one. Thus a = a,- and j3 = j3, will 
satisfy the goal. 

The idea is to apply this scaling for aU the tropical roots of tp{x) and each time, 
to compute n out of nd eigenvalues of the corresponding scaled matrix pencil, be- 
cause replacing P{X) by P{aiH) is expected to decrease the backward error for the 
eigenvalues of order a,, while possibly increasing the backward error for the other 
ones. 

More precisely, let ai < ai < . . . < denote the tropical roots of tp{x). Also 

let 




be the eigenvalues of P{pi) sorted by increasing modulus, computed by setting 
a = a, and j3 = t/:'(a,)^' and partitioned in d different groups. Now, we choose the 
ith group of n eigenvalues, multiply by a, and put in the list of computed eigenval- 



6 



Stephane Gaubert and Meisam Sharify 



ues. By applying this iteration for all / = 1 . . . c/, we will get the list of the eigenvalues 
of P{X). Taking into account this description, we arrive at Algorithm 1. It should 
be understood here that in the sequence /ii , . . ■,IJ.mi of eigenvalues above, only the 
eigenvalues of order a,- are hoped to be computed accurately. Indeed, in some ex- 
treme cases in which the tropical roots have very different orders of magnitude (as 
in the example shown in the introduction), the eigenvalues of order a,- turn out to be 
accurate whereas the groups of higher orders have some eigenvalues Inf or Nan. So, 
Algorithm 1 merges into a single picture several snapshots of the spectrum, each of 
them being accurate on a different part of the spectrum. 



Algorithm 1 Computing the eigenvalues using the tropical scaling 

INPUT: Matrix pencil P{X) 
OUTPUT: List of eigenvalues of P(A) 

1. Compute the corresponding tropical polynomial tp{x) 

2. Find the tropical roots of tp{x) 

3. For each tropical root such as a,- do 

3.1 Compute the tropical scaling based on a, 

3.2 Compute the eigenvalues using the QZ algorithm 
and sort them by increasing modulus 

3.3 Choose the ;th group of the eigenvalues 



To illustrate the algorithm, let P{X) = Aq +AiA +A2A^ be a quadratic polyno- 
mial matrix and let tp{X) = max(7o,7iA,72A^) be the tropical polynomial corre- 
sponding to this quadratic polynomial matrix. 

We refer to the tropical roots of tp{x) by a+ > a^. If a+ = which happens 

when 7^ < 7072 then, a = w ^ and j3 = tp{a)^^ = Jq^ . This case coincides with 



the scaling of [8] in which a* = 




When a+ ^ a , we will have two different scalings based on a+ = a =^ 
and two different /3 corresponding to the two tropical roots: 

^+ = tp(a+)-' = 5, j3-=t/,(a-)-' = l. 

rt yq 

To compute the eigenvalues of P(A) by using the first companion form linearization, 
we apply the scaling based on a+, which yields 

to compute the n biggest eigenvalues. We apply the scaling based on a^, which 
yields 
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to compute the n smallest eigenvalues. 

In general, let ai < ai < . . . < be the tropical roots of tp{x) counted with mul- 
tiplicities. To compute the ith biggest group of eigenvalues, we perform the scaling 
for a,, which yields the following linearization: 



V 



/l5af-'Ad_i ... i8a,Ai jSAqX 



V 








/ 



where j3 = t/?(a,) Doing the same for all the distinct tropical roots, we can com- 
pute all the eigenvalues. 

Remark 1. The interest of Algorithm 1 lies in the accuracy (since it allows us to solve 
instances in which the data have various order of magnitudes). Its inconvenient is to 
call several times (once for each distinct tropical eigenvalue, and so, at most d times) 
the QZ algorithm. However, we may partition the different tropical eigenvalues in 
groups consisting each of eigenvalues of the same order of magnitude, and then, the 
speed factor we would loose would be reduced to the number of different groups. 



5 Splitting of the eigenvalues in tropical groups 

In this section we state a simple theorem concerning the location of the eigenvalues 
of a quadratic polynomial matrix, showing that under a non degeneracy condition, 
the two tropical roots do provide the correct estimate of the modulus of the eigen- 
values. 

We shall need to compare spectra, which may be thought of as unordered sets, 
therefore, we define the following metric (eigenvalue variation), which appeared 
in [9]. We shall use the notation spec for the spectrum of a matrix or a pencil. 

Definition 1. Let Ai , . . . A„ and . . /i„ denote two sequences of complex numbers. 
The variation between A and /i is defined by 

v(A,^):=min{max|^„(,)-Ai|} , 

where S„ is the set of permutations of {1,2, . . . ,n}. If A,B G C"^", the eigenvalue 
variation of A and B is defined by v(A,Z?) := v(specA,specZ?). 

Recall that the quantity v(A,/i) can be computed in polynomial time by solving a 
bottleneck assignment problem. 
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We shall need the following theorem of Bathia, Eisner, and Krause [5]. 
Theorem 1 ([5]). LetA,B G C"^". Then v(A,fi) < 4 x 2-i/"(||A|| + |lB||)i-i/"||A - 

The following result shows that when the parameter 8 measuring the separation 
between the two tropical roots is sufficiently large, and when the matrices A2,Ai are 
well conditioned, then, there are precisely « eigenvalues of the order of the maxi- 
mal tropical root. By applying the same result to the reciprocal pencil, we deduce, 
under the same separation condition, that when Ai,Ao are well conditioned, there 
are precisely n eigenvalues of the order of the minimal tropical root. So, under such 
conditions, the tropical roots provide accurate a priori estimates of the order of the 
eigenvalues of the pencil. 

Theorem 2 (Tropical splitting of eigenvalues). Let f (A) = A^Ai + AAi +Ao 
where A; G C"^", and Ji := ||A,||, / = 0,1,2. Assume that the max-times polyno- 
mial p{X) = max(A272,A7i,7o) has two distinct tropical roots, a+ :— 71/72 and 
CC = Yo/Yi, and let 8 := a+/0£^. Assume that Ai is invertible. Let ^1, . . . de- 
note the eigenvalues of the pencil AA2 +Ai, and let us set ^„+\ — ■ ■ ■ = ^2n = 0. 
Then, 

v(specP,^)<^ , 

where 

C 4 X 2-'/2«(2 + 2condA2 + '-'/'"(condA2) , 

and 

a+(condAi)"' < |(^,| < a+condA2, 1 < / < « . (1) 

Proof. Let us make the scaling corresponding to the maximal tropical root 
hie 

Qin) = P+P{a+n) =A2Ai2+Ai^+Ao where 



Yi/YI' with j3+ = Y2/Y17 which amounts to considering the new polynomial matrix 



^2 = 72'^2> M^Yl^M, Ao^^Aq. 

Since A2 is invertible, A is an eigenvalue of the pencil P if and only if A = a+/Z 
where is an eigenvalue of the matrix: 



X = 



-A2 'Ai -At 'Ao 



/ 

Let jJ-iJ ~ 1 , . . . , 2« denote the eigenvalues of this matrix. Consider 

-A:7'Ai 
/ 
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Observe that ||Ai|| = 1 and ||Ao|| = JiYq/'^ ~ 1/^- Since the induced Euclidean 
norm 1 1 • 1 1 is an algebra norm, we get 

ll^ll < 11/11 + 114- 'Aill + IIA2 'Aoll < 1 + IIA2 'IIIIA2II + IIA2 'IIIIA2IIIIA0II 

= l+condA2(l + l/5) . 

Moreover, 

llrll < l+condA2 , = (condA2)/5. 

Using Theorem 1 , we deduce that 

v(specX,specy) <C/5'/^" . 

Since the family of eigenvalues of P coincide with a+(specX), and since the family 
of numbers coincides with a+(specy), the first part of the result is proved. 

If ^ is an eigenvalue of A2A +Ai, then, we can write ^ = o:+i^, where C, is 
an eigenvalue of A2/X +Ai. We deduce that 1^| < jjAj ' |1 |lAi || ~ condA2, which 
establishes the second inequality in (1). The first inequality is estabhshed along the 
same lines, by considering the reciprocal pencil of A2/X + A 1 . □ 

Remark 2. Theorem 2 is a typical, but special instance of a general class of results 
that we discuss in a further work. In particular, this theorem can be extended to 
matrix polynomials of an arbitrary degree, with a different proof technique. Indeed, 
the idea of the proof above works only for the two "extreme" groups of eigenval- 
ues, whereas in the degree d case, the eigenvalues are split in d groups (still under 
nondegeneracy conditions). Note also that the exponent in is suboptimal 

Remark 3. In [1, 2], the tropical eigenvalues are defined as follows. The permanent 
of an xn matrix B = (bij) with entries in Rmax is defined by 

perfi - max ^ bi„a^ . 

This is nothing than the value of the optimal assignment problem with weights {bij ) . 
The characteristic polynomial of a matrix C = (c,^) is defined as the map from Rmax 
to itself, 

xi^ Pcix) :=peT{C®xI) , 

where / is the max-plus identity matrix, with diagonal entries equal to and off- 
diagonal entries equal to —00. The sum C®xl is interpreted in the max-plus sense, 
so 

I max(c,i,x) if / = 



(C©x/),- 



The tropical eigenvalues are defined as the roots of the characteristic polynomial. 
The previous definition has an obvious generalization to the case of tropical matrix 
polynomials: if Co, ... ,Q/ are nxn matrices with entries in Rmax, the eigenvalues of 
the matrix polynomial C(x) := Co © Cix © • • • © C^/x^' are defined as the roots of the 
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polynomial function x i-^ per(C(x)). The roots of this function can be computed in 
polynomial time by 0{nd) calls to an optimal assignment solver (the case in which 
C{x) =Cq®xI was solved by Burkard and Butkovic [6]; the generalization to the 
degree d case was pointed out in [1]). When the matrices Aq, . . . are scalars, the 
logarithms of the tropical roots considered in the present paper are readily seen to 
coincide with the tropical eigenvalues of the pencil in which Q is the logarithm of 
the modulus of Aj., for <k <d. When these matrices are not scalars, in view of the 
asymptotic results of [1], the exponentials of the tropical eigenvalues are expected to 
provide more accurate estimates of the moduli of the complex roots. This alternative 
approach is the object of a further work, however, the comparative interest of the 
tropical roots considered here lies in their simplicity; they only depend on the norms 
of Aq, . . . ,Ad, and can be computed in linear time from these norms. They can also 
be used as a measure of ill-posedness of the problem (when the tropical roots have 
different orders of magnitude, the standard methods in general fail). 



6 Experimental Results 

6.1 Quadratic Polynomial Matrices 

Consider first P(A) = Aq + AiA + A2A^ and its linearization L ~ XX + Y. Let z 
be the eigenvector computed by applying the QZ algorithm to this linearization. 
Both ^1 = z(l : n) and ^2 = z{n + 1 ; 2n) are eigenvectors of P(A). We present our 
results for both of these eigenvectors; Tjj denotes the normwise backward error for 
the scaling of [8], and T], denotes the same quantity for the tropical scaling. 

Our first example coincides with Example 3 of [8] where IIA2II2 ~ 5.54 x 
10"^ IIA1II2 w 4.73 X 10^ IIA0II2 ~ 6.01 X 10-3 and A,- e £iOxW_ ^^^^ ^qO ran- 
domly generated pencils normalized to get the mentioned norms and we computed 
the average of the quantities mentioned in the following table for these pencils. Here 
we present the results for the 5 smallest eigenvalues, however for all the eigenvalues, 
the backward error computed by using the tropical scaling is of order 10^'^ which 
is the precision of the computation. The computations were carried out in SCILAB 
4.1.2. 



|A| 




'7(C2,A) 




'7.(C2,A) 




'7r(C2,A) 


2.98E-07 


l.OlE-06 


4.13E-08 


5.66E-09 


5.27E-10 


6.99E-16 


1.90E-16 


5.18E-07 


1.37E-07 


3.84E-08 


8.48E-10 


4.59E-10 


2.72E-16 


1.83E-16 


7.38E-07 


5.81E-08 


2.92E-08 


4.59E-10 


3.91E-10 


2.31E-16 


1.71E-16 


9.53E-07 


3.79E-08 


2.31E-08 


3.47E-10 


3.36E-10 


2.08E-16 


1.63E-16 


1.24E-06 


3.26E-08 


2.64E-08 


3.00E-10 


3.23E-10 


1.98E-16 


1.74E-16 



In the second example, we consider a matrix pencil with j|A2||2 ~ 10^^, ||Ai||2 ~ 
lO-', IIA0II2 ~ 10^ and A; e c40x40 Again, we use 100 randomly generated pencils 
with the mentioned norms and we compute the average of all the quantities pre- 
sented in the next table. We present the results for the 5 smallest eigenvalues. This 
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time, the computations shown are from MATLAB 7.3.0, actually, the results are in- 
sensitive to this choice, since the versions of MATLAB and SCILAB we used both 
rely on the QZ algorithm of Lapack hbrary (version 3.0). 



|A| 










J?7-(Ci,A) 




1.08E+01 
1.75E+01 
2.35E+01 
2.93E+01 
3.33E+01 


2.13E-13 
5.20E-14 
4.56E-14 
4.18E-14 
3.77E-14 


4.97E-15 
4.85E-15 
5.25E-15 
5.99E-15 
5.28E-15 


8.98E-12 
7.71E-13 
6.02E-13 
5.03E-13 
4.52E-13 


4.19E-13 
4.09E-13 
4.01E-13 
3.97E-13 
3.84E-13 


5.37E-15 
6.76E-16 
5.54E-16 
4.80E-16 
4.67E-16 


3.99E-16 
3.95E-16 
3.66E-16 
3.47E-16 
3.53E-16 



6.2 Polynomial Matrices of Degree d 



Consider now the polynomial matrix P(A) = Aq +AiX H l-A^/A'', and let L = 

X,X + Y be the first companion form linearization of this pencil. If z is an eigenvector 
for L then = z(l : «) is an eigenvector for P(A). In the following computations, 
we use to compute the normwise backward error of Matrix pencil, however this 
is possible to use any z{kn +1 : n{k+l)) for k = . . .d — 1. 

To illustrate our results, we apply the algorithm for 20 different randomly gener- 
ated matrix pencils and then compute the backward error for a specific eigenvalue of 
these matrix pencils. The 20 values x-axis, in Fig. 1 and 2, identify the random in- 
stance while the y-axis shows the logjo of backward error for a specific eigenvalue. 
Also we sort the eigenvalues in a decreasing order of their absolute value, so Ai is 
the maximum eigenvalue. 

We firstly consider the randomly generated matrix pencils of degree 5 where the 
order of magnitude of the Euclidean norm of A, is as follows: 



llAoll 


!|Ai|| 


11^211 


11^311 


IIA4II 


llAsll 


0(10-^) 


0(10-^) 


o(io-) 


0(10-') 


0(10-*) 


O(IO^) 



Fig. 1 shows the results for this case where the dotted line shows the backward er- 
ror without scaling and the solid line shows the backward error using the tropical 
scaling. We show the results for the minimum eigenvalue, the "central" 50* eigen- 
value and the maximum one from top to down. In particular, the picture at the top 
shows a dramatic improvement since the smallest of the eigenvalues is not com- 
puted accurately (backward error almost of order one) without the scaling, whereas 
for the biggest of the eigenvalues, the scaling typically improves the backward error 
by a factor 10. For the central eigenvalue, the improvement we get is intermediate. 
The second example concerns the randomly generated matrix pencil with degree 10 
while the order of the norm of the coefficient matrices are as follows: 



ll^oll 


ll^ill 


IIA2II 


11^311 


IIA4II 


IIA5II 


0(10-^) 


0(10-^) 


0(10-^') 


0(10-*) 
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In this example, the order of the norms differ from 10^^ to 10^ and the space di- 
mension of Ai is 8. Figure 2 shows the results for this case where the dotted line 
shows the backward error without scaling and the solid line shows the backward 
error using tropical scaling. Again we show the results for the minimum eigenvalue, 
the 40th eigenvalue and the maximum one from top to down. 





Fig. 1 Backward error for randomly generated Fig. 2 Backward error for randomly gener- 
matrix pencils with n = 20, d = 5. ated matrix pencils with n = 8, d = 10. 
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